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Abstract 



We study the electronic structure of a new type of Fibonacci superlattice 
based on Si <5-doped GaAs. Assuming that 5-doped layers are equally spaced, 
quasiperiodicity is introduced by selecting two different donor concentrations 
and arranging them according to the Fibonacci series along the growth di- 
rection. The one-electron potential due to 5-doping is obtained by means of 
the Thomas- Fermi approach. The resulting energy spectrum is then found 
by solving the corresponding effective-mass wave equation. We find that a 
self-similar spectrum can be seen in the band structure. Electronic transport 
properties of samples are also discussed and related to the degree of spatial 
localization of electronic envelope- functions. 

PACS number(s): 71.25.-s, 73.61.Ey, 71.45.Jp 
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I. INTRODUCTION 



Following the first fabrication of quasiper iodic semiconductor superlattices [ttlH, there 
has been a increasing interest in the study of one- dimensional systems describing quasiperi- 
odic structures. From the very beginning, most researchers have considered the Fibonacci 
sequence as a typical example of a quasiperiodic system HfD, and several characteristic prop- 
erties of Fibonacci systems have been reported during the last decade. Thus, it is now well 
established that Fibonacci lattices exhibit highly fragmented electron and phonon spectra 
with a hierarchy of splitting subbands displaying self-similar patterns ||, and their corre- 
sponding electronic density of states shows spiky features ||. This exotic electronic spectrum 
strongly influences electron propagation [^,|J and dc conductance through the system, even 
at finite temperature M. Furthermore, electronic wave functions are neither extended, in 
the Bloch sense, nor exponentially localised; instead, they are critical in Fibonacci lattices 

All these striking results, however, have been obtained for two particular kinds of 
model Hamiltonians, namely tight-binding models (either diagonal, off-diagonal or both) 
and Kronig-Penney models. Therefore, the question as to whether those features purported 
so far as characteristic of Fibonacci order are extensible to more realistic Fibonacci Hamil- 
tonians becomes very appealing from a theoretical point of view. On the other side, since it 
was realized that Fibonacci superlattices (FSLs) preserve its quasiperiodic order even if sub- 
stantial growth fluctuations in the sequential deposition of layers is allowed for 0, a consid- 
erable interest on the transport properties of quasiperiodic semiconductor heterostructures 
has emerged in the experimental field. However, FSLs studied up to date are quantum-well 
GaAs/AlAs based ones. This shortcoming does not seem necessary from an experimen- 
tal viewpoint, since computer controlled doping deposition techniques are currently able to 
construct a wide variety of superlattice potential profiles. In particular, epitaxial-growth 
techniques allow to prepare 5-doped semiconductor structures, in which a sheet of donors 
is localized within a few monolayers of the crystal. These impurity atoms supply electrons 
and give rise to strong confinement by the resulting one-electron potential. 

The main purpose of this letter is to show that distinctive features, previously obtained 
for simpler Fibonacci Hamiltonians, also appear in more complex and realistic systems, hence 
strongly suggesting that those features can be a universal fingerprint of one-dimensional Fi- 
bonacci systems. To this end, in the present work we propose a new type of FSL based on 
quasiperiodically Si 5-doped GaAs. We study the resulting electronic structure within the 
envelope-function and effective-mass approximation. The one-electron potential is found by 
solving the Thomas-Fermi equation. We find that the electronic spectrum is highly frag- 
mented and shows self-similar patterns which become very sensitive to the doping ratio 
concentration. The spectral type of our system is analysed by means of bandwidth-scaling 
techniques suggesting a underlying singular continuous character. Finally, transport prop- 
erties of the FSL at zero temperature are discussed in the context of the Landauer formula 
and related to the possible critical nature of the electronic states. 
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II. THE MODEL 



The system we study in this work is a FSL made of Si 5-doped GaAs. In general, a 
FSL of order N is generated from two basic units A and B by successive applications of 
the inflation rule A — > AB and B — ► A yielding a sequence of the form ABAABABA . . . 
This sequence comprises F/v_i elements A and i*V-2 elements B, Fi being the /th Fibonacci 
number given by the recurrent law Fi = F^% + i^_ 2 with the initial values F = F\ = 1. In 
the present model we take each unit (A or B) as a slab of GaAs of thickness a = 200 A with a 
Si 5-doped layer embedded in its center. The quasiperiodicity is then introduced by allowing 
the doping layers to take on two possible values of donor concentration, and arranging them 
according to the Fibonacci sequence. Each doping layer consists of a continuous positive 
slab of thickness d = 50 A with either or ionized donors per unit area. Hereafter 
we fix = 5 x 10 12 cm~ 2 , and allow to vary from 2.0 up to 6.5 x 10 12 cm~ 2 in 

order to study the influence of the doping ratio r = /N^ onto the resulting electronic 
spectrum. Finally, we assume that there exists a uniform p-type background doping with 
Na = 10 15 cm~ 3 acceptors per unit volume. This range of parameters corresponds to samples 
which could be easily grown by molecular beam epitaxy techniques now available. 

First we compute the one-electron potential Vj(x) in each basic unit (here j denotes 
A or B), where x indicates the spatial coordinate along the growth direction. This we 
accomplish by means of the Thomas-Fermi (TF) semiclassical model. It is known that the 
TF formulation is equivalent to the self-consistent formulation in a wide range of doping 
concentrations [JTTJ] , and it has been successfully applied to periodically Si 5-doped GaAs 



SLs [P^-[i^]. Hence we can confidently adopt this approach to calculate the one-electron 
potential due to the layer of ionized donors. Choosing the middle of the £-doped layer as 
the origin of coordinates, the nonlinear TF differential equation reads 

^-i^-v J („r- + ^.g + .).(i-) 

- SttA^^ + x^Q-z), j = A, B, (1) 

where ejr denotes the Fermi energy. The value of the Fermi energy cannot be computed 
self-consistently in the TF formulation, so it must be regarded as a parameter. Thus we 
fix the origin of energies so ep = henceafter. Distances and energies are scaled in units 
of the effective Bohr radius a* (= 100 A) and effective Rydberg Ry* (= 5.8 meV), respec- 
tively. The boundary conditions for this equation are those of the superlattice given by [T2 



Vj(0) = (a/2) = 0. The resulting potential becomes deeper in the gap on incresing donor 
concentration. 

We assume the validity of the effective-mass approximation, and we take isotropic and 
parabolic conduction band in the growth direction. This approach works well in some direct- 
gap compounds, as it is the case of GaAs. Once the potential Vj(x) in each basic unit is 
found from (|l|), the electron envelope-function and energy values can be obtained from the 
following one-dimensional Schrodinger equation 

- -^y + V(x)4j(x)=E^(x), (2) 
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where V(x) is the potential of the whole FSL. We assume that this potential at any point of 
the system is given by the one-electron potential Vj(x) we have just computed in each basic 
unit, where Va and Vb are arranged according to the Fibonacci series. Therefore, V(x) is 
simple a piecewise potential with F N regions. 



III. NUMERICAL ANALYSIS 

Since there is no analytical expression for the potential Vj(x), computations must rely 
on numerical procedures. To this end, we divide the each basic unit of the FSL in a grid of 
points {xk = kh}, where h = a/n is the integration step and n is the number of grid points 
in each basic unit. In our calculations we have taken n = 400, which is more than enough 
to obtain very accurate results. The discretized form of the Schrodinger equation (|2[) may 
be cast in the matrix form 

f ^(x k +i) \ _ ( ak -1 \ ( ip(x k ) \ = p f ip(xk) \ ( o\ 
{ *P(x k ) ) \ 1 ) ^(x fc _0 ) ~ n \ij{xk-i) ) ' [6) 

where we have defined a k = 2 + h 2 (V(xk) — E) for brevity. Iterating this equation one 
obtains 

( 1P(XN SL +1)\ _ p pi i>{x ) \ - T(N \ ( V»(»o) \ m 

T(Nsl) is the transfer matrix of the FSL and N$l = nF^ is the number of grid points in 
the whole structure. T(N SL ) is real and relates the envelope-function at both edges of the 
structure. Taking into account that T(k) = Pk T(k — 1) and T(0) = Pq we find the following 
recurrence relations involving only real parameters 

T u {k) = a k T 11 {k-l)-T 11 {k-2), 
T 12 (k) = a k T 12 (k-l)-T 12 (k-2), 
T 21 (k)=T n (k-l), 

T 22 (k) = T 12 (k-l), k = l,2---N SL . (5) 

These equations must be supplemented with the initial conditions Tij(-l) = 5y, Tn(0) = «o, 
^12(0) = —1, T 2 i(0) = 1 and T 22 (0) = 0. The resulting allowed energies can be found by 
imposing periodic boundary at both edges of the system with F^ layers arranged accord- 
ing to the Fibonacci sequence. Once we have obtained the transfer matrix for the whole 
FSL, the resulting energies for which the absolute value of the trace of the transfer-matrix 
corresponding to the whole FSL is smaller (larger) than 2 are allowed (forbbiden) [ |15|] . 



IV. RESULTS AND DISCUSSIONS 

In order to properly discuss the novel features arising from quasiperiodicity in our system, 
it is convenient to give a brief account of the electronic spectrum associated to this model 
when all 5 layers have the same donor concentration. In this case, which corresponds to a 
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usual periodic superlattice, resonant coupling between identical states of neighbouring layers 
leads to the formation of minibands of finite width [|T^ -1^|, related to extended (Bloch type) 



electronic states. The resulting miniband structure and the dispersion relation inside allowed 
minibands E(K), K being the momentum perpendicular to the layers, can be evaluated using 
the transfer-matrix approach as well. Since Bloch theorem must be satisfied in the periodic 



case, the dispersion relation is found to be (see Ref. |16j for details) 



cos(Ka) = -Tr [T(n)]. (6) 

Notice that the required time-reversal symmetry of the dispersion relation E(—K) = E(K) 
is conserved. As an illustration, we show in Figure |IJ the miniband structure for a particular 
realization corresponding to a donor concentration of 5 x 10 12 cm -2 . Note that the lower 
miniband is almost nondispersive since its width is rather small (about 0.4 Ry*). This 
suggests that the ground state of the Thomas-Fermi potential is only weakly coupled to its 
neighbouring wells. On the contrary, the second miniband is clearly dispersive, as seen in 
Fig. [I], and its width amounts 4.7 Ry*. The third miniband is even wider and it crosses the 
Fermi level. 

Now, we consider the most prominent features of the resulting electronic structure when 
quasiperiodicity is introduced. A schematic diagram of the FSL potential for F 5 = 8 is 
presented in Fig. ||]. From a mathematical point of view, one of the most characteristic 
properties of electronic spectra in Fibonacci systems is its highly fragmented, Cantor-like 
nature. We have confirmed this fragmentation in our FSL even when deviation from perfect 
periodicity is actually small, in other words, when the ratio r is close to unity. In fact, 
we have found that each miniband of the periodic SL, shown in Fig. [I], splits in several 
sub-minibands, that is, small gaps appear. The origin of these small inner gaps are directly 
related to the loss of long-range quantum coherence of the electrons, as the potential inside 
basic units, A and B, becomes different. Results corresponding to the fragmentation of the 
second miniband are shown in Fig. |asa function of the Fibonacci order N. We have mainly 
focused on this second miniband since the first one is almost nondispersive and investigation 
of the fragmentation process would require very tiny energy steps and then is rather time 
consuming. Only short approximants of the FSL are displayed in Fig. |3] since on increasing 
N the spectrum becomes so fragmented that it is difficult to observe minor features in the 
plot. However, we have carefully analyzed FSLs spectra up to order N = 12 (233 layers) 
and we have confirmed that the number of sub-minibands composing the whole spectrum is 
exactly F^, i.e. the number of basic units forming the superlattice. 

As we mentioned in the Introduction, another characteristic feature of Fibonacci systems 
is the self-similar pattern exhibited by their corresponding spectra. This self-similarity has 
been widely investigated within the tight-binding approximation while much less work has 
been devoted to quantum-well superlattices. Our results show that self-similar spectra are 
also obtained in Fibonacci Si 5-doped GaAs, as shown in Fig. (|. It is clear that the whole 
electronic spectrum for a short approximant (_F 4 = 5 in this case) is mapped onto a small 
portion of the spectrum of a higher approximant (F? = 21 in Fig. |4j). This is a consequence 



of how the FSL is constructed, based on a deterministic substitution sequence (17| . 

A third characteristic of Fibonacci systems concerns their spectral type, which results 
to be singular continuous. In order to estimate the spectral type associated to our model 
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Hamiltonian, we have calculated the normalized equivalent bandwidth S, defined as the 
ratio between the sum of all allowed sub-minibands and the width of the second miniband 
in the periodic A SL (4.7 Ry*). As can be expected from the Cantor-like nature of Fibonac- 
cian spectra, S vanishes as the system size grows. Furthermore, we have obtained that the 
normalized equivalent bandwidth decreases according to a power law of the form S = F^ 13 
with (3 ~ 0.1, as seen in Fig. [|. According with earlier works UTTf , such a behaviour is char- 



acteristic of a singular continuous spectrum for which all the envelope-functions are critical, 
i.e., regarding localization properties the functions are neither exponentially localized nor 
extended. 

The richness in structure displayed by the electronic spectrum should be reflected, to 
some extent, in its transport properties ||. We have evaluated the electrical resistance at 
zero temperature, p, using the well-known dimensionless single-channel Landauer formula 
P~8|1 : p = R(eF)/T(ep), where R and T denote the reflection and transmission coefficients, 



respectively. The transmission coefficient can be obtained in a straightforward manner in 



the transfer- matrix formalism [Kj. A typical example of the obtained results is shown 
in Fig. ^|(a) for Fn = 144. We observe that the Landauer resistance exhibits a highly 
fragmented structure displaying dramatic fluctuations under minor variations of the ratio 
r. The practical implications of such behaviour should be clear: One can select almost 
any desired value of the resistance of the sample by choosing properly the suitable ratio 
r. According to previous works 0, on increasing temperature interesting behaviours of the 
resistance can be expected since the detailed structure of the energy spectrum naturally 
determines the finer details of the resistance pattern at finite temperature. 

We have also investigated the spatial extent of the electron envelope-functions at the 
Fermi level. This we accomplished by means of the inverse participation ratio (IPR), as 
defined, for instance, in Ref. || 

IPR = Ek '^ )|4 2 , (7) 

(E fc M**)l 2 ) 

where the index k runs over grid points of the whole FSL. The IPR gives an estimation 
of the volume occupied by the electron envelope-function: The smaller the IPR, the more 
extended the electron state. Although a more detailed analysis of the envelope-functions is 
required to determine the exact nature of the wave function (multifractal analysis), for an 
understanding of the resistance behaviour the IPR is sufficient. Figure |6](b) shows the IPR 
for Fu = 144 as a function of the ratio r. Note that the minimum value is reached when 
this ratio becomes unity, i.e. in the periodic superlattice. The IPR also present dramatic 
fluctuations under minor variation of the donor concentration ■ A comparison of both 
Fig. |6|(a) and Fig. |6|(b) reveals that the IPR is larger (electronic states are more localized) 
whenever the resistance increases, indicating that the spatial extent of the envelope-functions 
controls the electrical transport of the sample. 



V. CONCLUSIONS 



We have proposed a new type of quasiperiodic (Fibonacci) superlattice heterostructure 
based on Si 5-doped GaAs. The corresponding electronic spectrum shows a highly frag- 
mented, self-similar nature resembling that found for simpler tight-binding models. The 
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spectral type of our model Hamiltonian, obtained from bandwidth-scaling considerations, 
indicates that it is singular continuous in the thermodynamical limit, in agreement with the 
current opinion of the mathematical community, supporting the conjecture that the spectral 



type for almost all substitution sequences should be singular continuous Hl9| . The spec- 
trum structure is very sensitive to the doping difference between the basic units, and this 
fact significantly affects the transport properties of the sample at zero temperature. This 
interesting result suggests the possibility of a certain degree of "engineering" of transport 
properties during superlattice growth by a proper selection of the corresponding doping se- 
quences. Finally, we are able to relate resistance fluctuations to the nature of the electronic 
wave function, through the inverse participation ratio. The obtained relationship suggests 
that the overall conductance of the superlattice is directly connected with the decay rate of 
the electron wave function along the sample. In this sense, it would be convenient to perform 
a complete analysis of the wave function espatial distribution by means of the multifractal 
formalism. Work in this direction is currently in progress and we expect to report on it 
elsewhere. 
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FIGURES 



FIG. 1. Miniband structure for periodically Si 5-doped GaAs with a donor concentration 
5 x 10 12 cm- 2 . 

FIG. 2. Schematic diagram of the one-electron potential for a Fibonacci Si 5-doped GaAs with 
= 5 x 10 12 cm" 2 and 7V"jf } = 4.5 x 10 12 cm" 2 with N = F 5 = 8 layers. 

FIG. 3. Allowed sub-minibands as a function of the Fibonacci order N, for a Fibonacci 
Si (5-doped GaAs with = 5 x 10 12 cm" 2 and = 4.5 x 10 12 cm- 2 . The number of 

sub-minibands is Fjv for each order N. 

FIG. 4. Self-similar spectrum of a Fibonacci Si 5-doped GaAs with = 5 x 10 12 cm -2 and 
= 4.5 x 10 12 cm -2 . Left plot shows the whole spectrum of a FSL of order A = 4 whereas 
right plot shows a detail of the spectrum of the FSL of order N = 7. 

FIG. 5. Log- log plot of the normalized equivalent bandwidth S as a function of the num- 
ber of (5-doped layers F/v for a Fibonacci Si 5-doped GaAs with = 5 x 10 12 cm~ 2 and 
= 4.5 x 10 12 cm- 2 . 

FIG. 6. (a) Landauer resistance at the Fermi energy and (b) IPR as a function of the ratio 
N^/N^ for a Fibonacci Si 5-doped GaAs with = 5 x 10 12 cm" 2 and F n = 144 layers. 
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